arXiv:1508.02855v2 [hep-ph] 5 Nov 2015 


TU-1002 

CERN-PH-TH-2015-189 


Higgs-bosons couplings to quarks and leptons 
in the supersymmetric Standard Model 
with a gauge singlet 


Andreas Crivellin^ and Youichi Yamada^ 

^ Theory Division, CERN, CH-1211 Geneva 23, Switzerland 
“^Department of Physics, Tohoku University, Sendai 980-8578, Japan 


Abstract 

The loop corrections to the couplings of Higgs bosons to quarks and charged leptons are 
calculated within supersymmetric versions of the Standard Model, extended by a gauge sin¬ 
glet. The effective couplings of the SU{2)l doublet and singlet Higgs bosons to quarks and 
leptons, induced by sfermion loops, are derived. Analytic expressions for the case of generic 
sfermion flavour mixing, including the complete resummation of all chirally-enhanced contri¬ 
butions are presented. These results are important in scenarios in which the mixing between 
singlet and doublet components of Higgs bosons is small, and the (pseudo) scalar compo¬ 
nent of the doublet is light. The calculated loop effects can have important consequences in 
flavour physics, especially for AF = 2 processes. 
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1 Introduction 


The supersymmetric model with a gauge singlet scalar in addition to the usual particle content 
of the Minimal Supersymmetric Standard Model (MSSM) is a well motivated model, and was 
already proposed in the early days m of studies of Supersymmetry (SUSY), as a solution 
to the fi problem [1]. It receives even more attention today in light of the recent Higgs boson 
discovery [5]. 

These models were for sometime even the focus of Higgs-boson studies in SUSY before the 
Higgs boson was actually discovered with the relatively low value of mass of ~ 126 GeV, which 
can be accounted for in the MSSM. The reason was the fact, noticed already quite early in 
Refs. [SHS], that the singlet coupling to the Higgs doublets allows to break the link of the 
Higgs qnartic couplings to gange conplings, typical of the MSSM. Supersymmetric models with 
additional gange singlets can therefore accommodate for a value of the SM-like Higgs mass much 
larger than that allowed in the MSSM, at least for small values of tan (3. Before the Higgs boson 
discovery, this fact reconciled SUSY aficionados with the possibility that the LHC would find a 
Higgs boson far heavier than that predicted by the MSSM (see for example Ref. m)- It also 
possess the capability of reducing the amount of parameter fine tuning needed in order to obtain 
the correct Higgs mass (see for example [nHi3]). Since then, it has become clear that a larger 
Higgs mass than that allowed in the MSSM is also possible for moderate to large large values 
of tan/3 mm due to the possible doublet-singlet mixing contributions to the physical Higgs 
masses. 

In addition to the original scale-invariant model, called the Next-to-Minimal Supersymmetric 
Standard Model (NMSSM), there exists other forms of the MSSM extended with a gauge singlet 
superfield. They differ by the singlet self interactions terms present in the superpotential. There 
is in particular the “Minimal Non-minimal Supersymmetric Standard Model” (MNSSM) |16lll7j 
or nMSSM, and the PQ-NMSSM [il (18112(1] . In all these cases, the only other particles to 
which the singlet superfield S couple are the Higgs doublet superfields Hu and H^. The models 
differ only in the way the Peccei-Quinn symmetry [/(1)pq [IlIISlET] of the superpotential is 
explicitly broken at the electroweak scale, i.e. either through a cubic term in S', or a linear 
one. The phenomenological studies of both models amounts to a sizable part of all beyond-SM 
analyses [SSHST]. 

Over the years, particular attention has been paid to the fact that the lightest pseudoscalar 
particle in these models, oi, can be quite light (see for example Refs. [281(31] . and references 
therein). Apart from small admixtures with the CP-odd neutral scalar components of Hu and 
Hd, the mass eigenstate ai is mainly the pseudoscalar component of the singlet snperfield S. 
The lightest CP-even neutral Higgs boson, hi, may be equally light, depending on the values 
of various parameters of the model. Potentially enormous consequence arise from this fact for 
Higgs studies at the LHC [29H32] . 

The presence of light particles in the spectrum clearly affects also the physics of flavoured 
mesons, modifying therefore the MSSM searches at the high-luminosity frontier. Indeed, the 
impact of very light Higgs bosons in K- and R-meson decays, has been the subject of intensive 
studies (see e.g. Refs. [23112411331(35] 1. Flavour physics will keep playing an important role in the 
exclusion or detection of snch light states throngh searches at the Belle H factory at SuperKEKB 
and with the LHCb program at CERN, in particular through the measurements of and 

mixing, and of the decay Bg —>• . Direct searches for ai at the LEP [361137] . 
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B factories [38VHT] . and the LHC [l2ll33] have already reduced the parameter space of these 
models. Needless to say, direct searches of additional Higgs states will be able to probe whether 
singlet states are part of the Higgs sector or not. 

It is therefore very important to know with a good precision the couplings of the Higgs bosons 
to quarks and leptons (of relevance for direct Higgs boson searches, flavour physics and also dark 
matter direct detection). Note that all Higgs mass eigenstates in this class of models, also those 
which are mainly singlet states, can couple to matter through their mixing with doublet Higgs 
states. Moreover, couplings of the singlet S to matter fermion, even though they are vanishing 
at tree-level, are generated at the one-loop level as was hrst pointed out in Refs. [441146] . and 
more recently in Ref. [35]. However, these effective singlet-fermion couplings have not attracted 
much attention. Although some existing codes for these models (such as NMSSMCALC [47] and 
SPheno [48] ) include these couplings partialljQ, no complete formula including resummation 
effects are implemented. It is true that the experimentally found value of the Higgs mass puts 
some emphasis on low values of tan j3, while these loop effects are maximized for the largest 
possible values of tan (3. Nevertheless, these effects cannot be neglected. They should be taken 
into account until the final embedding of the MSSM with a singlet is -hopefully- experimentally 
discovered. See, for example, Ref. [50] for the study of NMSSM with large tan/3. 

In this paper we plan to revisit the effective couplings of the singlet Higgs to leptons and 
quarks, induced by SUSY particle loopsJl To this end we work in the approximation of large 
tan /3 where the loop effects are phenomenologically relevant. This implies that a certain amount 
of hne tuning is needed to reproduce the measured value of the Higgs mass, like it is the case 
in the MSSM as well. With respect to the calculations in Refs. [44] . we drop the assumption 
of minimal flavour violation but rather include the effects of flavour-changing soft parameters. 
We also analytically re-sum all chirally-enhanced effects (as done in Ref. [51] for the MSSM) for 
threshold corrections to Yukawa couplings and to the CKM matrix. While also Ref. [46] worked 
in the MSSM with generic flavour structur^, we include the wino and bino contributions which 
are not included in Ref. [l6]. The results shown in this article are valid for all extension of the 
MSSM with a singlet, where the /r term is generated by the vacuum expectation value (VEV) 
of S, irrespectively of the singlet self-interaction terms allowed in the superpotentialEI 

This article is organized as follows: In Sec. 2, we review the basic properties of SUSY 
standard models with an additional gauge singlet supermultiplet. In section 3, the self energies 
of quarks and leptons generated by SUSY loops are calculated. Section 4 deals with the effective 
couplings of the singlet Higgs which are expressed in terms of these self energy. Section 5 shows 
our numerical results and hnally we conclude in Sec. 6. 

^In NMSSMTools [49] only the threshold corrections to the doublet Higgs Yukawa couplings but not the loop 
corrections to the singlet-quark couplings are implemented. 

^Note that the charged Higgs boson loops also give the tan/3 enhanced contribution in the NMSSM [3S]. We 
do not consider this contribution here, since it depends on the details of the Higgs potential. 

^Ref. |46| has presented numerical results only for the minimal flavour violation case. Note also that the 
calculation in Ref. [3S] has included the II — cl mixing. 

“^We assume a a CP-conserving Higgs potential. A recent study of the CP violation in the Higgs sector in 
NMSSM is seen in Refs. |52[I53| . 
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2 SUSY models with a gauge singlet 

In this section we review the basics of models obtained adding a singlet field to the MSSM 
particle content. In these models the superpotential is given by: 

w = + + f(S}, ( 1 ) 

where contains the usual Yukawa structure: 

lY^SSM ^ _ j^Cy^Q . ^ (2) 

with Q = {U,D)'^ and L = {N,E)'^, the quark and lepton SU{2)l doublet superfields; 
and E^, the SU{2)l singlet ones, and • symbolizing the antisymmetric SU{2) invariant product. 
The part of the superpotential denoted by f{S) contains the singlet self interactions that identify 
the specific model in this class. Before proceeding, we give also the soft SUSY-breaking scalar 
terms needed to specify these models: 

^soft = Yuk + Koft“ + ■ Hd + fsoitiS) , (3) 

with 

Koft“ = Q ■ Hu + Q ■ Ha + E^A^ L ■ Ha + H.c., (4) 

where the A-term couplings Af^ {F = U, D, L) are those often written as {E)AfjY^ in literature, 
and 


l/bmn ^ Q*mlQ + U^ml.U^*+ + L*mlL + E‘^m%.E^* 

H*uHu + H*aHa + m| S*S . (5) 

The part fsoit{S) contains the mass and self couplings of the singlet scalar. Throughout this 
article, the same symbol is used for the superfields Hu, Ha, S and their scalar component. In 
all the above expression flavour indices have been suppressed. We shall return to this issue later 
in this section. 

The gaugino mass terms: 

- ^ (Ms ^ + M2 WW + Ml Bb) , (6) 

together with the terms in Eq. ([3|) exhaust the list of soft SUSY-breaking terms in these models. 

As for parts of the potential describing the singlet self interaction, in the NMSSM, IU(5) is 
usually chosen to be 

IU(S) = ^«:53, (7) 

and Usoft(5') is therefore 

Usoft(5) = ^kA^S^ + H.c.. (8) 

A mass term +m'C<S'^ is at times added to Eq. m, and a corresponding one m'B qS"^ in Eq. ([8|) 
which break explicitly the scale invariance of the model. 

In the nMSSM, W{S) contains a linear term in S: 

W{S) = (+m'C5^), (9) 
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U{1)PQ 

U{1)r 

s 

-2 

2 

Hu 

I 

0 

Hd 

I 

0 

{u^Q) 

-I 

2 

(D^Q) 

-I 

2 

{EH) 

-I 

2 


Table 1: Peccei-Quinn and i?-symmetry charges of the various fields. 


where also an optional quadratic term in S, with a massive coupling m' can be added (see for 
example Ref. m), and 

VsoftiS) = m\C^S + {m'B(^S^) + H.c., (10) 

with and being massive parameters. 

No bilinear terms in are present in the above superpotential, nor in the soft scalar 

potential terms. It is assumed in these models that such terms are generated once S acquires a 
VEV Vs = (S). The effective fi and B parameters are then given by 


MefF 


Refr 


XVs, 

+ 


1 df*{S* 


dS* 


S=Vs 


( 11 ) 


respectively. 

Note that the mass term m? in Eq. Q may be generated spontaneously from non-renormalizable 
operators of the superpotential (in the Kim-Nilles mechanism [3]) or of the Kahler potential (in 
the Giudice-Masiero mechanism [SI]) involving the axion field, which acquires a VEV. There is 
in principle no particular constraint on this massive parameter m to be linked with the elec- 
troweak scale. A large tadpole term, however, tend to induce a very large value of vs- In such 
a case, an effective /ieg, compatible with electroweak-symmetry breaking, can be obtained only 
for a tiny value of A (which implies a very small singlet-doublet mixing). 

As is well known, the MSSM breaks explicitly the Peccei-Quinn (PQ) symmetry through 
the term. In the MSSM with a singlet it is possible to assign PQ charges to the various helds 
in such a way to have -|- A SHu ■ invariant under this symmetry (see first column in 

Tabled]). The explicit breaking of the PQ symmetry is thus shifted into f{S). 

It is clear that -I- A SHu ■ H^ also enjoys an R symmetry with charges for the various 

fields listed in the second column of Table dl We remind that the superpotential has i?-charge 
R{W) = 2. Scalar, fermionic, and auxiliary component of the same chiral superheld $ have 
different R-charges, with ii($ferm) = .^(‘^scai) ~ 1) and R(d'aux) = R(^scai) ~ 2. Moreover, 
the fermionic component \a of a gauge superheld 14, has R-charge R{Xa) = I- Thus, the R 
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symmetry must be broken in the process of breaking SUSY, generating the i?-violating mass 
terms for gauginos in Eq. (j6]) . 

Through the VEV of S, induced by the electroweak-symmetry breaking, both global U{1) 
symmetries are spontaneously broken. In particular, the breaking of the PQ symmetry results in 
an axion-like CP-odd scalar component of S', whose mass is proportional to n (or ^ depending on 
the specific model) and therefore vanishing in the limit k —)• 0 (or ^ 0). The limit k —^ 0 (or 
^ 0), in contrast, does not restore the R symmetry, because there are other terms in the 

Lagrangian that also break this symmetry. 

In addition to the PQ and R symmetric ones, there is also the MSSM limit, obtained for 
A, K —)■ 0 in the NMSSM (or A,^ —)■ 0 in the nMSSM) while keeping /rgff finite. In this limit, S 
tends to decouple from the model, giving rise to the superpotential and thus also to the scalar 
potential of the MSSM. However, the phenomenology may still be different from the one of the 
MSSM as the singlino (which is very weakly coupled for a finite but small A) can be very light, 
i.e. the LSP. Therefore, expectations for DM and collider searches may differ substantially from 
those of the MSSM. 

We close this overview of the models object of this paper with two subsections: one on the 
mixing of Higgs doublets and singlet, and one on the flavour basis we adopt in our investigation. 


2.1 Mixing of Higgs bosons 


The singlet Higgs boson and the neutral components of the doublet Higgs bosons mix with each 
other. In order to study this mixing we expand the singlet field S as 


5 = Us -h ^ {/is -k ios} , 
V2 

in analogy with and adopted in the MSSM 

~ Vd + {hd + idd} j 

Hu ~ W T 1 


( 12 ) 


(13) 


Here we have implicitly assumed that the global minimum of the Higgs scalar potential is realized 
by three real VEVs, u„, Vd and Vg- 

The presence of interaction terms among H^, and S in the scalar potential, such as 

ViH-s) = f^esX*{S + S*){vdH°a + VuH^) - XAxSivuH^d + ^dH^) + • • • , (14) 

induces mass-mixing terms among these components. The mass terms in the basis (|121 fT3]) are 
expressed as 

^mls? = - ^ h + a) , (15) 

with 

/ /id \ / \ 

h = I /i« I , a = I I . ( 16 ) 

\ hs / \ 
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The explicit expressions for the 3x3 mass matrices Mj and Mp in the NMSSM can be found, 
for example, in Ref. [25l[26]. 

The CP-even (CP-odd) mass eigenstates hi, /12 and (oi, 02 and unphysical Nambu- 
Goldstone boson G^) are given by 



where and are 3x3 orthogonal rotation matrices, such that 

{0^)'^MlO^ = diag(0,m2^,m2j. (18) 

Note that we choose our conventions in such a way that if the lightest CP-even (CP-odd) Higgs 
hi (ai) is mostly singlet-like we will have small angles in the mixing matrices. 

In the CP-odd Higgs sector, it is often convenient to use the “MSSM basis” as an intermediate 
step, to separate the Nambu-Goldstone mode G^ from physical states. The corresponding 
rotation is given by 

\ as J \ ai J 

with 

/ cos/3 sin/30\ / 1 0 

Op = I — sin /3 cos /3 0 j , = I 0 cos 9a sin 9a j , (20) 

\0 01/ \0— sin 9 a cos 9 a ) 

and ^ is the mass eigenstate in the MSSM. With just a rewriting of the above equations, the 
two pseudoscalars oi and 02 are then given by 



oi = sin 9aA^ + cos , 

02 = cos — sin 0yiOs . (21) 


Note that our angle 9a corresponds to the angle 9a — the existing literature. In the 

MSSM limit, when A —)• 0 with fixed /Tefr, we have 9a 0. The angle 9a, however, can be very 
small also in other cases. For example, in the Peccei-Quinn limit, in which oi is light, it is given 
by (see for example [26]) 


sin 9 a 


V sin (5 cos /3 

sin^ /3 cos^ /3 


( 22 ) 


which is greatly suppressed in the large tan /3 case. 


The rotation matrix for the CP-even Higgs states involve in general three different mixing 
angles. In the MSSM limit, it is 




— sin a cos a 0 
cos a sin a 0 
0 0 1 
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(23) 




with h 2 = h (the SM-like Higgs) and = H (the heavy CP-even Higgs), and the angle a 
defined as in the MSSM. In this same limit, if H, A, and are assumed to be heavy and 
nearly degenerate, the above matrix converges to: 




cos /3 sin (3 0 

sin /3 — cos P 0 

0 0 1 


(24) 


For reference, the physical state and the 
are the same as in the MSSM, 


Goldstone mode of the charged Higgs bosons 


( G^\^( cos/3 -sin/3 \ \ 

\ J V sin/3 cos/3 J \ J ' 


(25) 


2.2 Flavour parameters/violations 


We work in the super-CKM basis, in which the tree-level mass matrices, i. e. the Yukawa 
couplings, of the quarks and leptons in the superpotential are diagonal [5ll[55l[56]. Then the 
Yukawa coupling matrices in the superpotential have the following forms: 

{Yu)ij = (Yoh = {YL)ij = (26) 

where is the bare CKM matrix, arising from the misalignment between Yjj and Yd, which 
we have accommodated in the up-quark sector but one could have equally well shifted it into 
the down sector. The relation of the bare Yukawa couplings and to the physical 

fermion masses will be discussed in section 4. We also set (/ = u,d,i) to be real and 

positive. The doublet superfield Q is now defined as Q = D). In this basis, the tree- 

level couplings of neutral Higgs bosons to quarks are, indeed, flavour diagonal. Similarly, in the 
trilinear SUSY-breaking terms in Eq. (j4|), , A^, A^, have to be understood as 

(A'')., = A-iF/f, (A% = A/, , (Aq., = A'j, (27) 

where A^, and A^ are not necessarily diagonal and, in general, not even Hermitian. 

We now turn to the flavour violation of the sfermion (squarks and slepton) sector. The 
squark mass matrices in the super-CKM basis {qLi,QL 2 ,QL 3 , QRi,qR 2 ,qR 3 ) are parametrized as 


M} = 


2 2 
rrifTT Ta 


JLL 

(mj^R)^ m 


fRR/ 


(28) 


with Tn?fXY if ~ — [LL, RR, LR)) being 3x3 matrices in flavour space. The 

submatrices of Eq. (I28l) for down-type squarks are 


'^\ll 

= 

2 

~2 

^dRR 

= rriDc, 

i''^dLR)ij 

= -VdA\ 




Jij 


(29) 



f LFl 

Figure 1: The diagram shows our convention for the quark and lepton self energies SF . Here 
i and j are flavour indices and / = {u,d,i). 

For up-type squarks, we have 

mlLL = 

''^uRR = ! 

{mliRh = -VuA^* - (30) 

The mass matrix for the sleptons is obtained from by replacing [Q, D, A^) with (L, E, A^), 
repspectively. We drop 0{v'^^) terms in Eqs. [30]) in our approximation, as explained in 
Section 3. It should be noted that and related by SU (2) invariance. 


3 Quark and Lepton Self-Energies 


In this sections we calculate the complete set of chirality changing one-loop quark and charged- 
lepton self energies induced by SUSY particles (i.e. sfermions, gauginos and Higgsinos). 

We decompose the self-energies of quarks and leptons (see Fig. [T|) as 





Pr + 






Pl. 


(31) 


with incoming (SM) fermion fj and outgoing fermion /j. Here / = {u,d,i) denotes the fermion 
type and i and j are flavour indices. 

Since we know that the SUSY particles are much heavier than the SM fermions, it is possible 
to expand S^ (p) in powers of p/m-susY- For our purpose it is sufficient to evaluate the right- 
handed side of Eq. ([3T]l at p^ = 0, i.e. at leading order in p/msuSY- Eurthermore, since we are 
only interested in chirally-enhanced effects related to Higgs-fermion couplings we only need the 
chirality-changing part of the self-energies: 

= = (32) 

We further assume that the masses of the SUSY particles in the loops are sufficiently larger than 

f LR 

the VEVs of doublet Higgs bosons, i.e. {vd,Vu) <C msusY) and evaluate Sb to leading (first) 
order in v/M^usy- We refer to this approximation as the decoupling limit since the remaining 
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_jyia 


Figure 2: Quark self-energy with gluino and squark as virtual particles. 


terms of Sh do not vanish for (MsusY,/^eff) co. For calculating Sh to leading order in 
u/Msusy the S?7(2)-breaking elements of the SUSY mass matrices of the sfermions, neutralinos 
and charginos (such as the left-right mixing of sfermions or gaugino-higgsino mixing) are then 
not treated by the mixing matrices but rather by mass insertions involving or Vu- In this 
approach, all SUSY particles in the loops are the SU{2)l gauge eigenstates. We also need to 
drop the terms in the sfermion mass matrices, as is done in Eqs. [3n|l . in order to 

retain the non-decoupling terms only. 


The sfermion mass matrices still needs to be diagonalized due to possible flavour mixing in 
and Neglecting as explained above, the diagonalization is done as 


M) = 

= 


diag(m 


ft' 
Wf^ 0 
0 


2 2 2 2 2 > 

J2 JS Jl J2 J3 


(33) 


The 3x3 mixing matrices take into account the flavour mixing originating from the 

terms rnjn and respectively. Note that the relations m^L = m^L = fn-Lii = 1,2,3) and 

pj/nL _ y{o)y^dL fulfilled due to SU(2) invariance. 


For later convenience we introduce the abbreviations 

(34) 

where i, j, m = 1,2, 3. In Eq. (|34l) index m is not summed over. 


The self energies induced by sfermion loops resemble those in the MSSM by 
Note that the elements in Eqs. ()291130p . inserted into sfermion propagators in the loops, 

may generate chirality-enhanced effects with respect to the tree-level masses if they involve the 
large VEV Vu (tan /3-enhancement for down-quark/lepton self-energies) or a trilinear ^■I-term 
(2l{^./(Y/MsusY)-enhancement). 

Below we list the relevant contribution (as calculated in Ref. m)- 

We start from the gluino-squark contributions (see Fig. [2]). In our approximation, this 
contribution is proportional to and given by 


LR 


2as 


3 


^ imiLR)jk Kki^ CQ[mj, 

m,n=l 


j,k=l 


Hm 


mjR 

dR 


(35) 
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Figure 3: Fermion (quark and lepton) self-energy with sfermions and neutralinos as virtual 
particles. 





Figure 4: Fermion self-energies with sfermions and charginos as virtual particles. 


Here Cq is the standard three-point function 


at vanishing momenta: 


= (70(0,0, 0; m3) 


2 _ 2 \ 


3/ 


d^k- 


TT^ J — m\ + ie\ [A:^ — m^ -|- ie\ \kP‘ — m| -|- ze] 

_ m\ m| + m| m| ln(m|/m|) -|- m| m^ In(m 3 /m 3 ) 

{m\ — — m'^){rri^ — m^) 


(36) 


For the neutralino-sfermion contributions to lepton and quark self-energies shown in Fig. [3] 
we get 


^ex° LR _ 


1 


IGvr^ 

3 


+ E 


m=l 


dlMi {mlLR)jk ) Cq 

^ j,k=l m,n=l 
1 


2 2 2 
, rrifL , mm 


Mw sin P (^glM2fleS Cq (^\M2\‘^ , \gesf , 


+ 


V^92 

5lA7l/ieffC'o (|A7ip , l/ieffp )) 

g2^MK sin/3Mi/reff yC(o)aI^R Cq (jM^f , , m 




(37) 
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^dx° LR 


1 


lOvr^ 
3 






j,k=l m,n=l 


+ 


+ 


^ sin/? ( 5 iM 2 /xeff Co 

1=1 LV y? 


\ldes\ ,rngL 

Hm 


al 


A^l/leffCo ^|Afi|^ , l/ieffl^ , 


ym 


+ ^al aeS Y^f Co (I Ml I ^ , I /ieff I ^ rn% 


dR 


(38) 


Finally the chargino-sfermion contributions to lepton and down-quark self-energy (see Fig.[2|) 
are given by 


s^dx^ LR 
^fi 


^£x^ LR 
Yfi 


ydi(O) 

IGvr^ 


Meff 


<5.3 ^ qf* Fj™ („4„)53 a: 


L u RR 
^n33 


m,rx=l 


m=l 


V2g2sml3MwM2 ^ 

^yb(O) 


IGvr^ 


-Aiefrfi'2 sm/3MwM2 ^ A^^^jCofm?^ , |/iefip , IM 2 I 


m=l 




(39) 


In Eqs. (I38l[3^ . 0{vd) terms of the gaugino-higgsino mixing are neglected since they do not lead 
to chirally-enhanced contributions and cause unnecessary complication due to their ultraviolet 
divergences. In addition, in the higgsino-squark-squark contribution of Eq. (I39p we have further 
neglected small up-type Yukawa couplings of the first two generations and multiple flavour- 
changes, instead of the full form 


LR 


IGvr^ 


'Meff 


E 

/'d'd=i 


ET 




jj 


T/(0) \dLL 



(40) 

By using this approximation, we can hnd an analytic resummation formula taking into account 
all chirally-enhanced corrections m- 


In contrast to down-type quarks, the up-type quark self-energies cannot be enhanced 

by tan/3. Nevertheless, an enhancement by Af^/Y^Msusv is possible for the gluino and bino 
diagrams. These contributions are given as 






2 a. 


37r 


' E E i^lLR)jk Kki^ 

j,k,j',f'=l m,n=l 


16^ E .(41) 


m,n=l 
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We then denote the sum of all contributions as 

LR 

s^u LR 

f LR 

As to be discussed later, the flavour off-diagonal pieces of Sgenerate the flavour-changing 
couplings of the neutral Higgs bosons to quarks and leptons. In the case of the “minimal 
flavour violation” (MFV) with flavour-diagonal sfermion mass matrices, only charged higgsino 
contribution in Eq. ()39l) cause flavour mixing originating from the CKM matrix. 


jndg LR LR LR 




"fi 


LR lr 




fi 


Y^glR Y'^X' 


°LR 


( 42 ) 


4 Renormalization and threshold corrections 


As already stated, in the fermion self-energies (I35H39P the Yukawa couplings of fermions 

should be understood as the running ones of the superpotential in the singlet-extended SUSY SM. 
These couplings can be calculated from the physical masses of the fermions by properly taking 
into account the SUSY threshold corrections. As is well known, the chirally-enhanced corrections 
to the down-type quark masses may become numerically significant at large tan /3 |551I59I[64] and 
must be resummed to all orders. In addition, the off-diagonal self-energies cause rotation 
of the fermion mass eigenstates in the flavour space, which also generate difference between the 
bare CKM matrix of the superpotential and the physical one V. 

In this section, we review the procedure of obtaining yf(o) y(0) f^Qj^ g]yj running 
masses m/. and CKM matrix Vij, including the resummation of the chirally-enhanced correc¬ 
tions, following the results of Ref. |51j . 

The running mass rUg. of the quark qi extracted from experiment using the SM prescription, 
is given by 

m,, = ^, {q = u,d). (43) 

Q LR 

Here is the flavour-diagonal piece of the self-energy calculated in the previous section. 

Note that all terms in Eq. (1431) have to be evaluated at the same renormalization scale, i.e. the 
SUSY scale. 

In the down-type quark sector, is decomposed into the part which is proportional to 

a Yukawa coupling and the one which does not involve a Yukawa coupling, as 


■^dLR 


\^dLR 






Y 


di{0) 


(44) 


This decomposition is possible if we restrict ourselves to the decoupling limit where we have 
terms proportional to one power of Y'^ho) most as can been see from Eq. (13511 . (I38p . and (|39p . 
The second term of Eq. ()44p gives chirally-enhanced corrections to the relation between the 
quark masses and the Yukawa couplings of the superpotential, i.e. modifies this relation via 
a chirally-enhanced threshold correction [SgEHMSS- One automatically resums all chirally- 

®For a 2-loop analysis of the threshold corrections to the relation between the Yukawa couplings and the quark 
masses see Ref. |65 |. 
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(45) 


enhanced corrections by inserting Eq. 


into Eq. (H3]) and solving for 


Vd{l + tan/3ef) 


The corresponding expressing for leptons follows trivially by replacing d with i. In contrast, 
no such resummation is necessary for E“do) -v^^here only the contribution from ^d-terms can be 
significant: 

_ yuLR 

y“i(o) =_ ^ _ ('40) 


f LR 

We now turn to the flavour-changing part of Sj.- which modifies the relation between the 
physical CKM matrix and the CKM matrix of the superpotential. In order to simplify the 
notation it is useful to define the quantity 


ai = 


^fLR 


(47) 


max{mf.,mf^}' 

f LR 

for i ^ j. The elements Sj- contribute to the fermion mass matrices and therefore require an 
additional rotation with respect to the super-CKM basis to obtain the physical mass eigenstates 
of the fermions 

^fLiR)^^fLiR)^fLiR)^ ( 48 ) 

To leading order in small ratios of the quark masses mf^/mf^ <C 1, ^ then reads |66I168] 


/ 


= 


a- 


12 


G 


/ \ 


13 


—a 




f* 

12 


a. 


23 


J* 


I <^13 (Tiq CJ. 




12 "^23 


—G, 


/* 

23 


1 


(49) 




G 


The corresponding expressions for U-^ ^ are obtained from the ones for U-^ ^ by the replacement 

I ^ Gf\ 

Applying the rotations in Eq. (I49jl to the undjiW^ vertex renormalizes the CKM matrix. 
The bare CKM matrix in Eq. (j26l) can be calculated in terms of the physical CKM matrix 
V as 

yW = U^LyjjdLt _ (50) 

However, we have to take into account that in Eq. (1501) depends on through the 
chargino loop contribution to In general, an iteration procedure is necessary to calculate 

using Eq. (I50p . Nevertheless, in our approximation, we find a closed form of We first 
decompose G'j- as 

i=3 


— 

^fi - 


I d y(0)*TA(0) 

^f3 ^ ^FC'^Sf 


33 ! 


G 


fi 


i=l,2 


(51) 


SO that Gj- does not depend on (off-diagonal) CKM elements and 
ydsfO) 3 


^FC — 


-1 


E Aif,{mlLR)hK§3^Co(M\rn^ • (52) 

777 , , \ Hm ^71J 


m,n=l 
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In terms of the generalized Wolfenstein parametrization defined in the appendix of Ref. m , we 
find 


with 


1--h ^ Wim 


1 /( 0 ) = 


-Vh 

^12^23 - ^13 


1 - e 


die 

FC 


1 - 


1 ^ 12 ! 


V 12 

2 


^23 


l-£' 


dir: 

FC 


^^13 


\ 


-tvim 


1 _ 

^23 

1 _ p^ 

1 tp(^ 


(53) 


fl 2 — ^^12 + ~ ^12 ) 


^^23 — '^’23 + <723 — (T23 


Vis — Vis + <^13 — <^13 + <^ 12^23 + [^12 ~ ^12j^23 ~ ^ 12<^23 ’ 


Vjm = V12 Im 


^12 + ^^12 


— Im 


'Zdie 

^ 12^12 


(54) 


and 


/1- 


V = 


i-|^;i2r/2 


^>12 


Vl 3 


—V 


12 


^12^^23 - ^ 


1 -\V12\^/2 V 23 


13 


—U. 


23 


(55) 


f LR 

We are now in a position to relate the self-energies Sh to corrections to the Higgs-fermion 
vertices and compute the effective Higges couplings as a function of SUSY breaking terms. 


5 Effective Higgs couplings to quarks and leptons 

In this section we compute the effective Higgs-fermion couplings in the MSSM with an additional 
gauge singlet superfield. For this purpose we first determine the couplings to the Higgs doublets 
{Hu and Hd) and the singlet 5, and go afterwards to the physical basis with diagonal Higgs 
mass matrices. 


5.1 Couplings to doublet Higgs bosons 


The calculation of the effective couplings of the Higgs doublets Hu and H^ is the same as in the 
MSSM. We first decompose as 


rp(d,e) LR 
yiU LR 


Tp{d,t) I Tpl{d,t) 

El- ’Vd + E-) 'v.. 
E^^Vu + E'^^Vd, . 


(56) 

(57) 


The first terms on the right-handed side of Eq. (j56p is the holomorphic part, generated by the 
loop correction to coupling fufjRH^ which exists already at the tree-level and is only induced 
by H-terms in our approximation. In contrast, the second non-holomorphic term of Eq. (I56p 
comes from the loop-generated effective fiifjRHu coupling involving the effective /Xgff term as 
shown in Fig. [5j The term in Eq. (j56p proportional to Vu is always accompanied by a factor 
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Figure 5: Decomposition of the down-type quark self-energy into its holomorphic part Efj and 

its non-holomorphic part This decomposition is possible in the decoupling limit, i.e. if 

f LR 

Sh is evaluated at leading order in u/msusY- For charged leptons one simply has to replace 
d with i. 

Vs due to the PQ symmetry. For Eq. (f57)l . similar discussion holds by exchanging {Hd,Vd) and 
{Hu^Vy). Note that Ejj and E[ - are general functions of Vs- 

The loop contributions to the effective couplings to doublet Higgs are then obtained by 
replacing Vd and Vy in Eq. (l56|) (Eq. (IFT)! ) hy {Hd) and Hy {H*), respectively. 

5.2 Effective singlet-fermion couplings 

We are now in a position to calculate the loop-induced couplings of the singlet Higgs S to 
quarks and charged leptons resumming all chirally-enhanced corrections. Note that also these 
loop-corrections are directly related to quark and lepton self-energy generated once the singlet 
acquires its VEV Vg giving rise to figg. Therefore, the effective couplings in our approach can 
also be calculated in terms of the chirally-enhanced quark self-energies as done before for the 
doublet Higgs couplings. For the couplings to the singlet we only consider down-type quarks and 
charged leptons since the contributions to up-type quark couplings are not chirally enhanced. 
We hrst derive the couplings in the super-CKM basis and relate them to the effective coupling 
in the physical basis with diagonal quark and lepton mass matrices. The effective Lagrangian 
governing the interactions of quarks and leptons with the singlet is given by 

LeS = '^qf^i^hsqfLQiR + i^q^gi^asqfLQiR + h.c. + {q^i). (58) 


5.2.1 Gluino and bino induced couplings 

For the effective quark-singlet coupling induced by a gluino shown in Fig. [6] it is sufficient to 
replace in Eq. (p5]l by A in order to get the corresponding effective coupling to the singlet. 
The same is true for the pure bino contribution to the quark or lepton self-energy. The expression 
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{df)L 



Figure 6: Feynman diagram showing the gluino contribution to the effective singlet-down-quark 
couplings in the decoupling limit (i.e. in leading order of v/msusy)- 


for the effective gluino induced coupling then reads: 


pg/is LR _ pgos LR _ f '^'*^11 j^fdg 

dfdi dfdi /* 


(59) 


Here the hat refers to the fact that the couplings are given in the super-CKM basis. The 
expressions for the bino contribution to effective quark or lepton couplings is simply obtained 
by exchanging the corresponding sub- and super-scripts. 


5.2.2 Chargino and neutralino induced couplings 


Concerning the effective couplings induced by loop-diagrams with chargino and the neutralinos 
the situation is more involved since //eff can also appear from the diagrams like the ones shown 
in Fig. [71 This effect is contained to all orders within the higgsino propagator emerging from 
the Dyson series 


11 11 1^11 1*1 1 

- ,3 _ I,,2 

AAAAAAAAAA ^1 ^eff 


(60) 


Each /ieff = Aus arises from a coupling to the singlet. In the corresponding self-energy in the 
decoupling limit (Eq. pH]) and Eq. pp) i we pick out only the part of the propagator proportional 
to /ieff (or /r*jj) meaning that we necessarily have an odd number of couplings to the singlet whose 
contributions to the Dyson series is: 


1111*11 MefT 


(61) 


In all except one of these couplings the singlet is replaced by its VEV Vg. If the coupling is to 
S (and not S'*) the Dyson series gives 


1 1 1 ^ 

J\j J\, /v /t /V 


k‘^\S 


XS 




-f 


(62) 
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Figure 7: Example Feynman diagrams showing the chargino contribution to the effective singlet- 
down-quark couplings in the decoupling limit (i.e. in leading order of v/msusy)- 


where the factor 2 takes into 
other hand, we have instead 


account the possible permutations. If the coupling is to S*, on the 




(63) 


Thus we can write for the chargino induced quark-singlet coupling in the following way: 


LR _ 

1 

dfdi ~ 

"71 

\X^aaLR _ 

1 

dfdi ~ 

"71 




Meff 


fi 


Vu T-^/dx^ 




d 


(64) 


heff 


Of course also the lepton-singlet coupling induced by charginos can be written in the same way 
and the corresponding formula for the part of the neutralino self-energy which contains gaugino- 
higgsino mixing is straightforward. Note that the formulas (I64p are also valid for gluino and 
bino contributions, where the second term vanishes. 

We denote the sum of all contributions as: 

t^clsLR _ . T^x^asLR . Y\gasLR j^hsLR _ j^X^hsLR . j^^hsLR . j^ghsL/R 

'-dfdi -'^dfdi +^dfdi +^dfdi > '-dfdi -'^dfdi +^dfdi +^dfdi ■ 

and for leptons d is simply replaced by i. 


5.3 Higgs couplings in the physical basis 

Until now, we calculated the effective couplings of the Higgs doublets and the singlet to matter 
fermion. For this we worked in the interaction eigenstate for the Higgs sector and in the super- 
CKM basis with diagonal Yukawa couplings for the quarks and leptons. Both for quarks/leptons 
and for the Higgses, this is not the physical basis with diagonal mass matrices. Therefore, 
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additional rotations (see Eq. (jl7l) and Eq. (I49p i are required. Note that since the masses of the 
quarks and leptons entirely originate from the doublets Hu and their relation to the Yukawa 
couplings is the same as in the MSSM. 

After integrating out the heavy SUSY particles, interactions of the doublet Higgs bosons 
with quarks are described by the Lagrangian 


/•sCKM _ /oa 

^doublet ~ L 






Tpfd Tja 

^fi^u 


d-i 


R 


+ Q^f 


fL 


Y-6fi + E-AebaHt* - 




UiR + h.c. , 


( 66 ) 


with and E'^ determined by Eqs. ()56p and ()57p . After EW symmetry breaking the quark 
mass matrices are given by: 


m% = Vd [Yf6f, + E%>l + VuE% , (67) 

+ VdE% . (68) 


We use these relations in order to eliminate the dependence on Y^ and in Eq. ()66p . In 
addition, we go to the physical basis with diagonal quark mass matrices 


Vl 


R 


Lz-k q j tQ Jri 

jf ^jk^ki =mq^Ofi. 


(69) 


In the physical basis with diagonal quark mass matrices, the doublet-Higgs interactions with 
quarks are given by 


■^doublet — dj R 

- UfL 
+ "dfLVfj 
+ dfLV*f 

Here we defined the quantitj@ 


Vd 


^ + f/i diu 


m. 


■Ui 


Sf,-eYcot(3]H^u^ + eYH: 


ttO 


fi -'-'d 


Vd 


5ji — (cot /3 + tan /3) e^j 


m. 


Ui 


6ji — (tan /3 + cot /3) e“j 


diR 
UiR 

Ej- diR 
HEviR + h.c. 


(70) 


= {Ud^^E'dud^)fi. (71) 

which gives rise to the chirally-enhanced corrections to the Higgs couplings in the physical basis. 
In the Lagrangian 


Leff = Tq}q’:^HlqfLqiR + iTf^q';''alqfLqiR 

+T^^d^^H+UfLdiR + T^^u!:''H-dsLUiR + {h.c.) , (72) 

®We have dropped most of the parts contributing to and therefore £“, in Eq. (HU, since their contributions 
to the singlet Higgs couplings are suppressed by cot fi and irrelevant for our study. See, for example, Refs. |51l69l70j 
for the contributions of e“ to the effective couplings of 77^ to up-type quarks. 
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with = {hd,hu) and this leads to the following effective Higgs couplings to 

down-type quarks 


p/irf LR 

dfdi 


LR 

dfdi 


LR 

^ dfdi 


^sin/3(!^(i^,-eJ,tan/3) 


(73) 


while the couplings of the Nambu-Goldstone boson ^df<k^ ~ (1/v^) cos f3{mdi/vd)Sfi receive 
no corrections from e'^. 

For leptons the neutral Higgs vertices follow trivially from the ones for down-type quarks. 

The quark held rotations in Eq. (1691) also lead to a redehnition of the quark-singlet couplings 
from the super-CKM basis as 


■y^hgi^CLs) LR _ jjdL*‘-r^^s{s^s^ djRjjdR 

^dfdi -^rf^df,d,, ^i'i ^ 

■p/ls(as)-L-R _ -^'-Rrr^R 

^ifii ' 


(74) 

(75) 


Finally, we obtain the following couplings to the Higgs mass eigenstates: 


/ t^L2 LR 

^ dfdi 

LR 

dfdi 

■p/ll LR 

\ ^dfdi 


{o^f 


( ^ 
r 

V r 


hd LR 

dfdi 

h-n LR, 

dfdi 

hs LR 

dfdi 


( r 


G°LR 

dfdi 

T^ci2 LR 

dfdi 

LR 
df di 



/ -pG°LR 
'^dfdi 

YA° LR 

^ dfdi 

ptts LR 

\ ^ dfdi 


(76) 


The analogous results for charged leptons follow by replacing d with i. 

At this point we can already make a rough estimation of the effective singlet couplings to 
down-type quarks (and charged leptons as well). Comparing Eq. (fM)l and Eq. (173]) . it 

is seen that the singlet couplings are suppressed hy v/vg compared to the loop-induced couplings 
of the doublets {hd, A°). As a result, the singlet couplings decouple with (Msusy, A^eff) oo 
and hxed A, while the effective doublet Higgs boson couplings remain hnite in the same limit. 
Furthermore, the scaling with tan fd is different: while the loop-induced parts of the doublet 
couplings scale as tan^ ,0, the singlet couplings scale as tan/3. Nevertheless, for the Higgs mass 
eigenstates which are almost pure singlet, as is the case for example for pseudo-axions in Peccei- 
Quinn symmetric limit (I22j) . may give dominant contribution to the effective couplings to 

down-type quarks. 


6 Numerical Results 

In this section we study the numerical significance of loop-induced effective couplings of the 
singlet Higgs to down-type quarks and charged leptons. We first illustrate the generic size 
of these couplings and then discuss the potential effect of the singlet-induced contributions to 
flavour physics. 
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In the following we will quantify the size of the effective couplings of singlet Higgs bosons 
to down-type quarks (charged leptons) For the couplings to quarks we 

include QCD running and evaluate them at the squark mass scale. As input values for SM 
parameters (a^, nii,, ...) we use the current PDG values m- For the SUSY sector parameters, 
we use the following values: Mi = M 2 = 1 TeV, rrig = 3 TeV, A = 1. The chirality-conserving 
submatrix of the sfermion mass matrix in Eq. (j28p is assumed to be diagonal and flavour- 
independent, rriq = 2 TeV and = 1 TeV for squarks and sleptons, respectively. In the 
discussion of the flavour-changing couplings, however, we introduce the mixing between the third 
and second generations in the left-left or right-right sector. The chirality-conserving submatrix 
is then expressed as 


2 2 

^dLL = rriq 


10 0 \ 

0 1 , 
0 1 J 


(77) 


while and (/ = it, d, €] are obtained by replacing by and 62 ^^, respectively 

(and rriq for sleptons). The A-parameters in the chirality-changing submatrices are set to 


(A“)ij = diag(0,0,A*), {A^)ij = diag(0,0,A'’), {A^)ij = diag(0,0,A^). (78) 

We use A^ = A"^ = 1 TeV. Other parameters are specified in the figure captions. As stated 
before, we consider the case where CP violation is negligible and set all parameters to real. 


6.1 Flavour-conserving couplings 

We first show the flavour-conserving effective couplings of the singlet Higgs bosons to the down- 
type quarks and charged leptons, for the case of minimal flavour violation (MFV) — 

^ 23 ^ = 0). Since these couplings are purely loop-induced, they are much smaller than the 
couplings to the doublet Higgs, which exists already at the tree-level. This is especially the case 
for charged leptons where the loop corrections cannot involve the strong interaction. 

As shown in Eq. (j64h . there are two types of the contributions to the singlet couplings: the 
one directly proportional to the non-holomorphic part of the self energy E'^ and the other 
involving the derivative with respect to /iefr- For CP-conserving case, the latter only contribute 
to the couplings of the CP-even component hg, causing the different behavior between 
and 

Fig. El shows the loop-induced singlet couplings to tau leptons. As expected, these couplings 
are much smaller than the Higgs-tau coupling in SM, mT-/(\/2u) = 0.0072. 

For down-type quarks, the behavior of the effective couplings is similar to that for leptons, 
but two additional numerically important parameters enter (beside that evidently squark masses 
and slepton masses enter): A^ and rUg. Furthermore, the threshold correction to the relation 
between and mf, is much larger than that for r lepton (because is involved) which leads 

to an asymmetric behavior of with respect to the sign of //eff- We show the dependence 

of the effective singlet coupling to bottom quarks qj^ the (A*, /Xeg) plane in Fig. [9l 

Again, they are smaller than the coupling of the SM Higgs, mb/{V^v){at 2 TeV)~ 0.0095. 
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Figure 8: Left: Strength of the loop-induced coupling of tau leptons to the singlet Higgs for 
tan/3 = 50 and A = 1 as a function of the (effective) /x parameter. Blue (dark gray): hs, red 
(light gray): a*. Right: Same as the left plot for tan/3 = 30. 



Figure 9: Left: Strength of the loop-induced coupling of bottom quarks to the singlet Higgs 
bosons, to hs (left plot) and to (right plot), for tan/3 = 50 and A = 1 in the /r — plane. 
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Figure 10; Left: Strength of the loop-induced coupling of tau and muon to the singlet Higgs 
induced by (^ 23 ^^ (and normalized to 623 "^) tan/I = 50 and A = 1 as a function of the 
(effective) fi parameter. Blue (dark gray): hg, red (light gray); a^. 

Right: Same as left plot for tan /3 = 30. 

6.2 Flavour-changing couplings 

The flavour-changing couplings of fermions to the singlet Higgs and to the doublet Higgs are 
both loop-induced. As shown in Section 5, the singlet couplings are suppressed by the factor 
~ (v/vg) cot (3 compare with the corresponding doublet couplings. Here we focus on the flavour¬ 
changing couplings between the second and third generations, induced either by the flavour 
mixing of left-handed sfermions <523 or by the CKM matrices in chargino-sfermion loops. 

6.2.1 Leptons 

The plots of Fig. [TOl show the behavior of effective coupling of the singlet Higgs induced 

by 523 "^ ■ contribution is to a very good approximation proportional to unless ^ 23 "^ 

very large. The effect of 62 ^^ to jg always suppressed by the ratio m^/rrir compared 

to the ^ 2 ^^ contribution, and therefore in most scenarios subleading. This is also the case for 
the contribution of to since it does not involve the W loops. 


6.2.2 Quarks 

Let us consider first the case of the MFV { 62 ^ = 0). In this case only the quark-squark-chargino 
vertex induces flavour-violation. In Fig. [11] we show the size of the effective singlet couplings 
induced via chargino loops. 

In the presence of non-minimal sources of flavour-violation in the squark sector, additional 
contributions are induced. In case of a non-vanishing element Yi^s>o-s)LR -g ggnerated. 

Since in this case no interference with the MFV contribution occurs, their sizes are proportional 
to ^ 23 ^^. In the presence of 62 ^^, in contrast, there is interference with the MFV contributions 
for as seen in Fig. [12] 
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Figure 11: Size of the effective singlet couplings to bottom and strange quarks for MFV in the 
plane, for tan /3 = 50 and A = 1. 
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Figure 12: Size of the effective singlet couplings to bottom and strange quarks in the 
plane for // = 2 TeV, tan j3 = 50 and A = 1. 
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Figure 13: Contributions to the mass difference in the Bg system generated by chargino box 
diagrams (blue) and by the singlet exchange (red) for tan/3 = 50 (left plot) and tan/3 = 30 
(right plot) for A = 1. We assume the singlet-like has a mass of 7 GeV while hg and doublet 
Higgs bosons are much heavier. We set = 3 GeV. 

6.3 Effects in flavour-changing processes 

In general, the flavour-changing neutral current processes which are most sensitive to the ex¬ 
change of neutral scalars which couple to quarks proportionally to their mass are Bg ^ —)• 
and Bg^d — Bg^d mixing. For large tan/3, these processes receive potentially large contributions 
from the double-penguin diagrams mediated by neutral Higgs bosons [72H71]. However, as dis¬ 
cussed previously, the effective (flavour-conserving) couplings of the singlet Higgs scalars to 
charged leptons are much smaller than the ones of the doublet Higgs scalars. This means that 
the possible effect of the singlet exchange in Bg d —> must be much smaller than the effect 

of the doublet Higgs exchange. For Bg^d — Bg^d mixing, in contrast, two loop-induced flavour¬ 
changing Higgs-quark couplings are involved in both the doublet and singlet exchanges. As a 
result, singlet contributions scales compared to the doublet contribution as tan^ /3) 

{niA is the typical mass of the heavier doublet Higgs bosons), which can be non-negligible for 
very light Qg. Furthermore, in the MFV case the doublet Higgs contribution to Bg^d — Bg^d 
mixing tends to chancel between the GP-even and GP-odd states having similar masses IzuEa], 
while such a suppression is not efficient for the singlet states where hg and Ug can have very 
different masses. 

We illustrate the relative importance of the singlet exchange in Fig. [13] for the MFV case and 
nia^ = 7 GeV. The singlet exchange contribution is evaluated by using the effective quark-singlet 
couplings at the renormalization scale rua^. The meson form factors are adopted from Flavour 
Lattice Averaging Group (FLAG) [75]. We see that the singlet Higgs contribution can be much 
larger than the chargino-squark box contribution. 


7 Conclusion 

In the framework of the MSSM extended by a gauge-singlet supermultiplet (e.g the NMSSM), 
we have studied the loop-induced effective couplings of the singlet Higgs bosons to quarks and 
charged leptons. Allowing for the most general flavour structure of the SUSY breaking terms, 
we have derived analytic formula for the couplings of the singlet to quarks and leptons. Keeping 
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the leading term in the expansion of f/MsusY we have performed the complete resummation 
of all chirally-enhanced effects, making our results valid in the phenomenologically large tan/3 
region. 

We found that loop-induced singlet couplings to down-type quarks and charged leptons 
are enhanced by one power of tan/3, while the flavour-changing couplings of down-quarks to 
heavy doublet Higgs bosons are enhanced by tan^ /3. In addition, the loop-induced singlet- 
quark couplings vanish in the decoupling limit v <C {Msusy , fJ-es) while the doublet couplings 
remain finite. Nevertheless, the loop-induced singlet couplings can be the dominant part of the 
couplings of the lightest states oi and/or hi to quarks and charged leptons if they are to a 
good approximation singlet like. Furthermore, these couplings can also be phenomenologically 
important if such ai is very light as it give enhanced contributions to flavour observables. 

While an analysis of the impact of the effective couplings of the singlet Higgs in realistic 
scenarios, for example in the NMSSM, is beyond the scope of this article, we have pointed out 
that sizable effects of singlet Higgs in flavour observables are still possible. In our numerical 
analysis we have examined the generic size of the effective singlet quark (lepton) couplings, both 
flavour-conserving and flavour-changing ones, and considered their impact on processes where 
the singlet contribution is particularly relevant: for very low singlet masses, we have shown that 
it can be the dominant beyond-SM contribution to the Bg — Bg mixing due to the enhancement 
by light mass of ag. 
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